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ABSTRACT 

We present here a theoretical model to account for the stellar IMF as a result of 
the composite behaviour of the gas and dust distribution functions. Each of these 
has previously been modelled and the models tested against observations. The model 
presented here implies a relation between the characteristic size of the dust grains 
and the characteristic final mass of the stars formed within the clouds containing the 
grains, folded with the relation between the mass of a gas cloud and the characteristic 
mass of the stars formed within it. The physical effects of dust grain size are due to 
equilibrium relations between the efhciency of grains in cooling the clouds, which is a 
falling function of grain size, and the efficiency of grains in catalyzing the production 
of molecular hydrogen, which is a rising function of grain size. We show that folding 
in the effects of grain distribution can yield a reasonable quantitative account of the 
IMF, while gas cloud mass function alone cannot do so. 
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1 INTRODUCTION 

When we observe the mass distribution function of stars in 
a galaxy we obtain the present day mass function (PDMF) 
which is not the mass function of stars at their birth (the 
initial mass function, IMF) but the distribution resulting 
from the accumulation of stars since the first star formation 
occurred in the Galaxy. In this process, whether the star 
formation is continuous or sporadic, the most massive stars 
with their short lifetimes are more quickly eliminated, and 
the masses of all the stars are reduced secularly by mass 
loss, while the stellar population may be augmented by the 
arrival of stars from outside the Galaxy in accretion events 
with dwarf galaxies. The local IMF is a basic function in the 
study of the chemical evolution of our Galaxy, as it is the 
measure of the number of stars formed in a given mass in- 



terval in the solar neighborhood (|Scalolll978l . Il986l : iTinslevI 
1 19801 ). The origin of the IMF is a fundamental problem in the 
whole of astrophysics because it determines the photomet- 
ric properties of galaxies, and the dynamica l and chemical 
evolution of their interstellar media. ISalpeteJ (|l955l ) was the 
first to note that there appears to be a power-law relation- 
ship between the number of stars observed in the Galactic 
field and their masses, and it was his work which gave rise 
to our concept of the IMF. 
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In order to try to account for the form of the IMF one 
needs to know the mechanism of the process of star for- 
mation. We understand that this process begins with the 
formation of giant molecular clouds within which, given the 
right c onditions of t emperature and density (the Jeans con- 
ditions IJeand (|l902l )) local gravitational collapse gives rise 
to stars. It is also well understood that the presence of dust 
grains as catalysts is an essential condition for the conver- 
sion of atomic to molecular hydrogen, so that they are nec- 
essary for the presence of the giant molecular clouds. It thus 
appears reasonable that in order to investigate the factors 
which explain the stellar IMF one needs to take into ac- 
count not only the gas cloud mass distribution function, but 
also the distribution function of the dust grain sizes. The 
relationship between the gas and star formation is obvious 
since stars form from the gas clouds, but the connection be- 
tween the dust and star formation is less direct. As we have 
just noted, the presence of dust, even in small quantities, 
makes a major differe nce to the con version rate of atomic 
to molecular hvdrogen ISchavj (|2004l ). but dust also affects 
the heating cycle in the ISM, which affects the tendency to 
form the cold dense clouds needed for star formation. 

Arguments based on gravitational instability and on ob- 
servations of molecular gas reveal that the low-temperature 
(T ~ 10 K) high-density (n ^ 40 cm"^) cores in giant molec- 
ular clouds (GMG's) are the natural sites for stars to form. 
Although individual GMC's are well resolved in the Milky 
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Way and in our local group (|Bolatto et alj l|2008h . and ref- 
erences therein), it is from CO observations of nearby spi- 
rals which show most clearly that star formation occu rs in 
regions dominated by molecular gas ( | Wong fc Blita i2002l : 
iKennicutt et al.ll2007i : iBigiel et al.ll20m 



The present article, whose aim is to bring out the con- 
nection between cloud mass distributions, dust grain size 
distributions, and the stellar mass distribution, is organized 
as follows: in Section 2 we present the version we will be 
using of the stellar IMF derived from published observa- 
tions by a number of authors, in Section 3 we develop three 
models in which the observed mass functions for the Gi- 
ant Molecular Clouds in the Galaxy will be used, in con- 
junction with the interstellar dust grain size distribution, 
to predict the IMF: the zero model (which we will combine 
with the three models of gas clouds to obtain the IMF) is 
our numerically derived dust grain s ize di stribution func- 
tion taken from ICasuso fc BeckmanI (|2010l ) (CBIO) which 
fits well the observed carbonate and silicate distributions, 
and takes into account both the production and the destruc- 
tion of grains, the latter by grain-grain collisions in the ISM. 
The first model uses the numerically der ived molecular cloud 
mass distribution function, taken from ICasuso fc BeckmanI 
l|2007t l (CB07) in which clouds may coagulate to form big- 
ger clouds or may be disrupted in collisions, depending on 
their masses, temperatures, densities and relative velocities, 
while a fraction condense to form stars. The second model 
uses an analytically deri ved molecular cloud mass distribu- 
tion function taken from lCasuso Sz BeckmanI (|2002i ) (CB02) 
where we solve the stochastic differential Langevin equation 
in a context of existing barriers due to box-effect. The third 
model uses a fully analytic approximation based on the as- 
sumption of balance between the thermal emission by dust 
and dust heating due to collisions with gas molecules, and 
is based on the CB02 approximation, in Section 4 we com- 
pare the predictions of the models with the observationally 
derived IMF taken form several authors, and in Section 5 
we present our conclusions. 



2 THE OBSERVATIONALLY DERIVED IMF'S 

As the observational references required to test our models 
we have used three observational studies of the local stellar 
mass distribution at birth (i. e . the IMF). The first refer- 
ence is lMaciel fc Rocha-Pintol (|l998l ) who derive two possi- 
ble IMF's, one of which assumes a constant star formation 
rate (SFR) and the other assumes a variable SFR derived 
from observations of late-type stellar populations on the 
main sequence. Both of these versions of the IMF used the 
presen t day stellar mass function (PDMF) bv lKroupa et al.l 
l|l993ll . We have plotted these in Fig. 1, where we can see 
the difference in the slopes of the two curves due to the well 
understood trend of stars to disappear from the PDMF at 
a greater rate the higher the initial mass, and the different 
effect this produces using different histories for the SFR. 
The change in slope required to produce the IMF from the 
PDMF is greater with the variable SFR, because this in- 
clude s a relatively recent observed peak in the local Galacti c 
SFR (|Maciel fc Rocha-Pintolll998l : lRocha-Pinto et al.ll2000D . 

We have use d two oth e r refer ences for the Galactic IMF, 
one obtained by iKroupal l|200lh which takes into account 



Salpeter (1955) 
Kroupa (2001) 
PDMF Kroupa (1993) 

Maciel & Rocha-Pinto (1998, variable SFR) 
Parravano et al. (2010) 
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Figure 1. Observationally based stellar mass functions for the 
solar neighbourhood, compared with the classical Salpeter IMF. 
X(log M) is the number of stars per logarithmic mass 
interval. 



constraints from local star count data, an d IScald (|l998l ) com- 
pilation of MF power law indices for young clusters and OB 
associations, and which shows different slopes for different 
mass ranges, and an other, the semi-empirica l Galactic IMF 
obtained recently bv jParravano et"ah 

I (|201ll ). both of them 

also shown in Fig. 1. 



3 THE THEORETICAL MODELS, FOLDING 
THE GAS AND DUST DISTRIBUTION 
FUNCTIONS 

3.1 Semi-analytical treatments 

3.1.1 The numerically derived dust grain size distribution 
function (GSDF) 

Knowing that ISM dust grains are generally formed in the 
outer atmospheres of stars in the late stages of their evolu- 
tion (mainly red giants, asymptotic giant branch stars, and 
SNe) and that grain lifetimes can be estimated as less than 
100 Myr, we have taken an observational version of the local 
stellar birthrate to calculate the radius distribution of grains 
surviving today. To do so, we computed the mass distribu- 
tions of stars present over this period assuming a standard 
IMF, the observed stellar birthrate function, and stellar life- 
times as a function of their masses, and folded in the initial 
distribution of dust grain sizes at each epoch, according to 
a simple prescription of grains produced in different stellar 
mass ranges. The differences between the size distribution 
of grains produced in the older, low mass stellar popula- 
tion and those produced in younger, high mass stars include 
the tendency of the smaller grains to be preferentially de- 
stroyed by shocks in the post-SN environment. Under these 
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assumptions, the major peak in the grain size distribution 
comes from dust produced in the younger more massive 
stars, and the three peaks in the interstellar dust grain size 
distribution function (GSDF) of carbonaceous grains corre- 
sp ond to the three maxima in the local SFR, as detailed 
in iRocha-Pinto et alj ()200Cl[ ). This correlation between the 
three main star formation events in the Gala:xy and the 
three peaks observed for the sizes of grains appears clear. 
For silicates, the observed GSDF shows only one peak, close 
to the biggest peak in the distribution observed for carbona- 
ceous grains. This increases the plausibility of the scenario in 
which the peak for carbonaceous grains at the smallest radii 
is associated with those grains produced during an interval 
of some 10 Myr near the late-stages stars of lower masses 
(near 1 M0 arising from the earliest major peak in the SFR) 
because in that case the velocities of grains in the expand- 
ing shells are lower than in stars with higher masses so that 
the time interval for collisions among grains is greater so the 
probability of shattering is higher and the mean sizes of the 
grains are lower. For the second peak, we have grains com- 
ing from stars of intermediate age (intermediate mass com- 
ing from the second peak in the SFR) where the velocities 
of grains are intermediate and so some of the larger grains 
can escape from the high density shells where the shatter- 
ing can break them into smaller grains. Finally the third 
peak, with the largest grains, can come from SNe where the 
expansion velocities are so high that the largest grains can 
escape from the shells into the ISM where the densities are so 
low that the collisions occur relatively infrequently. The ob- 
served difference in the slope of the distribution function for 
carbonaceous (close to 0) and silicate gr ains (close to 1 ) may 
be due to the mechanism proposed by IVidali et all (|2005l ') 
for suprathermal grains, whereby the mislaignment with re- 
spect to the interstellar magnetic field, may be stronger for 
carbonaceous dust than for silicate dust. This would lead to 
greater coUisional effects for the carbonaceous grains, chang- 
ing the initial slope to that observed. The global features of 
the GSDF are rather well reproduced by the model in which 
we use only the dust sizes yielded by the stellar production 
process and the accompanying modifications due to shatter- 
ing in the local environment of the late stages of the stars, 
together with a stellar birthrate function from observations. 
However, a closer fit to the carbonaceous grain distribution 
at the low radius end is obtained if we also include the mod- 
eling of longer term grain-grain collisions processes in the 
ISM (CBIO). In Fig. 2 we can see the complete model (full 
line) and the comparison with data. 



3.1.2 The numerically derived molecular cloud gas mass 
distribution function (MCGMF). 

As one can see in CB07 we parameterize the distribution 
mass function of gas clouds at a given epoch as a function of 
their initial distributions of density, temperature, velocity, 
and mass. The position and velocity distributions are three- 
dimensional. We assume different sets of initial distribution 
functions, (Gaussian, flat, power-law) for temperature, den- 
sities, masses, and velocities of an initial sample of ~1000 
gas clouds, either pressure bound or gravitationally bound, 
evolving within a confinement volume specified as a cubic 
box, and arranged initially as knots in a cubic grid. Because 
our best fit to data is found for Gaussian distributions, and 
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Figure 2. Model of the GSDF produced by combining the stellar 
grain production curve with the effects of size modification due 
to collisions in the circumstellar environment. Observational data 
sets are shown for comparison. 



for spheroidal gas clouds we have taken that case for the 
present work. Each cloud has an effective volume that de- 
pends on its mass and density via . After a given evo- 
lution time, with a given initial set of velocities, we consider 
that two of the clouds have collided if the distance between 
their centers is less than or equal to the sum of their radii r, 
which are calculated using r = (f^)^''^- To take particular 
values for the main physical parameters (density, tempera- 
ture, velocity, and mass) we take the expression: 



1 



cr(27r)i/2 



2l/2„ 



(1) 



where f{x) is the probability distribution function of a vari- 
able X centred on a value Xo and with a dispersion a. We 
have inverted these functions using the formula: 



X ^Xo + AX(-2a'/o<?[/i(a;)])cos[27r/2(j:)] 



(2) 



where /i and /2 are two random values of the distribu- 
tion, and AX is the range of values permitted for the 
variable X. To take into account the phenomenon of gas 
cloud collapse to form stars, we assume that those gas 
clouds that have attained the Jeans mass disappear and 
we count them as new stars (however, if these stars have 
masses big enough that their lifetimes are less than the time 
during which the model runs, these stars also disappear). 
We have taken for the Jeans mass the classical expression 
Mj = i7r^/^p~^/^«jG~^/^, where p is the variable density 
of the gas cloud, Vs is the sound speed, and G is the grav- 
itational constant. If t — ti> ^^'^,9° then a star with mass m 
formed at time ti is not computed at the current time t. The 
input free parameters of the model are: (1) the width and 
centre of the Gaussians (normalized to 1) for temperature, 
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velocity component (one for each of the three coordinates), 
mass, and density; (2) the ranges of values permitted for 
temperature, velocity component, mass, and density; (3) the 
critical values adopted for temperature, velocity and density, 
and of course the number of time steps considered for the 
time evolution. At each time step, if the temperatures of the 
two colliding gas clouds are greater than the adopted critical 
temperature, the relative velocity of the clouds is less than 
the critical velocity, and the densities of the two clouds are 
less than the critical density, i.e., the collision is completely 
inelastic, then we assume that one cloud disappears, while 
the other cloud change to a new cloud with a mass equal to 
the sum of the masses of the initial two clouds (the clouds 
merge). If all the previous conditions occur, but the relative 
velocity between clouds is greater than the critical velocity, 
then we assume that both clouds disappear by diffusion into 
the diffuse ISM. If the temperatures of the two clouds are 
less than the critical temperature, and the two densities are 
greater than the critical density, then we differentiate two 
cases: (1) when the relative velocity is less than the critical 
velocity, then star formation occurs and the two clouds dis- 
appear; and (2) when the relative velocity is greater than 
the critical velocity, we assume that each of the two clouds 
breaks into two subclouds, for simplicity each subcloud hav- 
ing half of the mass of the parent cloud. If the temperature 
of one cloud is less than the critical temperature and the 
temperature of the other is greater than the critical temper- 
ature, and the density of the first cloud is greater than the 
critical density while the other density is less than the crit- 
ical density, then the first cloud remains unchanged, while 
the second cloud breaks into two subclouds each with a mass 
half the mass of the parent cloud. All the kinematics are 
computed subject to momentum conservation in each kind 
of collision. So for the first case, i.e., a completely inelastic 
collision (merger) , the output velocity for each of the three 

dimensions is computed as v'^ = *^m'i+M2''' ' where the su- 
perscripts 1 and 2 indicate each of the two colliding clouds, 
and the subscript i indicates each direction in space (x, y, 
z). For the case in which each cloud breaks into two sub- 
clouds, we assume that for the four resultant subclumps, 
those assigned 1 and 3 have the same x, y velocity com- 

i i 1 • / 2M2V2 i , (Mi—M2)vi i , 

ponents computed via Vi = j^PTm^ + hh+M^ ' ' where 
Vj^i indicates the velocity component i of the initial cloud 
j, and the z-direction will be such that v'^^ ., — —v'^ .,, with 
v'i,z = ^'z- Similarly for the other two output subclumps 2 

and 4, usmg Vi = j^^^ m.+m^ ' ' '"2,^ = -«4,., 

with U2,z = v'z. For the case when one cloud remains un- 
altered (called 1), while the other breaks into two pieces 
(called 2 and 3), the output velocity of cloud 1 changes to the 
-, while the other two components 



2M2V2^i (Mi — M2)vi 



M1+M2 



M1+M2 



(Mi-M2)u2,, 



2 and 3 have the output velocities v'^ = M1+M2 — 

and 1)2,2 = ~^3,z- To simulate the fact that the clouds lie on 
gravitationally bound orbits, we impose box shaped bound- 
ary conditions, such that the sign of the relevant velocity 
component changes when the cloud reaches the edge of the 
box. One can see in Fig. 3 the result compared with data. 




logM(Msun) 



Figure 3. Model of the Galactic gas cloud mass distribution 
function, taken f rom CB07, (sol i d line ) compared with observa- 
tional data from lSolomon et al. dotted line. The form of 
this function is not similar to that of the stellar IMFs shown in 
Fig. 1, but does play a role in determining the IMF as explained 
in Section 3 of this article. 



3.1.3 The analytically derived molecular cloud gas mass 
distribution function. 

The spaces between the gas clumps inside a cloud are large, 
i.e., they show a low filling factor of 0.1 for clumps within gi- 
ant molecular clouds (GMCs), implying a distance between 
clumps of ~6 pc while the size of each clump is typically of 
1 pc. A similar structure is observed for the clouds in the 
galactic disk as a whole: characteristic distances between 
clouds are of 5 x 10^ pc while the size of a cloud is of order 
10 pc. These observational data, together with the observed 
relative velocities: 5 km s~^ between clumps and 50 km s~^ 
between GMCs in the disk, permit us to make the approxi- 
mation that the general distribution of gas in the ISM may 
be driven by random collisions leading to a Brownian (Gaus- 
sian) distribution function of displacements, at least after 
a sufficiently long period of time. In fact, we can see how 
the clumps of gas clouds have number distributions against 
size which in log-log plots appear as perfect Gaussian distri- 
butions (see CB02). Translated to linear units this implies 
Planckian distributions divided by size, implying an origin 
based on a system formally equivalent to standing waves in a 
box, although in fact somewhat more complex. As a first ap- 
proximation we consider the stochastic differential Langevin 
equation for the time evolution of a sample of clouds under 
Brownian motion and also suffering the effects of a magnetic 
field B. Assuming that B is uniform, time independent, and 
perpendicular to the velocity it of clouds, for simplicity, one 
has: 



(3) 
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where q and m denote the charge and mass of the cloud re- 
spectively. In this equation the influence of the surrounding 
medium on the cloud can be split into three terms: firstly, 
a systematic term —Xu representing a friction experienced 
by the cloud, due to its movement through the other clouds; 
second, a fluctuating term A{t) which is the characteristic of 
stochastic Brownian motion; and third, the simplified influ- 
ence of an assumed uniform and constant magnetic field B. 
We have taken a magnetic force proportional to the velocity, 
because the magnetic force per unit volume at a distance R 
from the galactic centre is , and the magnetic flux in a 

galactic disk is puocB^ . The frictional term —Xu is assumed 
to be governed by Stokes' law: the frictional force decel- 
erating a spherical cloud of radius b and mass m is given 
by ^'^^^ , where rj denotes the coefficient of viscosity of the 
surrounding fluid. Ifence one has A — 22!*1, xhe fluctuating 
term A{t) has been restricted for simplicity to be indepen- 
dent of u and with variations extremely rapid compared to 
the variations in u. To solve equation (3) first we write it as 



(see CB02). Selecting the expression for a single wall and 
taking xo = and tio = we have: 



with 



P = X + 



qB Qnbr] -\- qB 



Now the solution to equation (4) when t^fi ^ is: 



W{x,t;XQ,UQ) 



(4) 
(5) 

(6) 



a kind of diffusion coeficient, and 



with D = jTTT^ 
W{x,t;xo,uo) is the probability of displacements of length 
X — Xo in any given direction at time t starting from ini- 
tial values Xo and uo- To obtain the same result from an- 
other physical point of view, we can consider the clouds 
as aggregates of microclouds. Then, the probability that 
a microcloud arrives at a cloud surface after a path of 
length xi (assumed approximately the same as the distance 

between clouds) at time t, is P{xi,t) — 2t(Dtl)^/^ ^ 
IChandra8ekha^(|l943l '). and so, assuming that the size of mi- 
croclouds is approximately constant, the distribution of sizes 
X for the final clouds should be W{x,t) oc f^ P(x\,t)dx\, 
i.e., a Gaussian distribution similar to equation (6). To make 
the change of variable from displacements x to masses m, we 
assume that the distribution of displacements of clouds be- 
fore collision with other clouds of similar size, is comparable 
to that of the sizes of the clouds. This assumption, taken to- 
gether with the observational relation between the sizes and 
masses: xocm" with 0.27 ^5 a^C0.41 leads to a distribution: 

dN ^js^ 

- — ocm e " (7) 
dm 

However, because the regions where the clouds move (mainly 
the spiral arms) are regions where the motion is limited spa- 
tially, we can treat the motion formally as if were confined 
within a box with walls which, although not in fact fixed, 
may be treated as fixed to a first approximation. And this 
situation can be formulated by a distribution function of 
displacements which is the addition of several Gaussians: 
the main Gaussian (centred on zero), and the other Gaus- 
sians centred on the walls which reflect the clouds whose 
displacements are distributed following the main Gaussian 



W{x,t) 



+ 6" 



(8) 



where L/2 is the size of the region, and the maximum of 
W{x,t) is at Xo = L/2. We have neglected losses of clouds 
through the barrier (i.e. out of the confined region). Since 
xoam" , one can assume LolMq, and then one has the final 
distribution: 



dN 
dm 



(9) 



3.1.^ Folding the gas and dust distribution functions. 

We propose here that there is a monotonic dependence of the 
mass scale of stars forming within a given gas cloud, and the 
scale of the dust grain size within the cloud. This hypothesis 
has a priori plausibility because it is based on the known 
requirement of dust grain surfaces to catalyze the forma- 
tion of H2 molecules from an initial cloud of HI. The fractal 
nature of the dust grains implies that the larger are the 
grains the greater is their effective catalytic surface to vol- 
ume ratio, so that larger grains tend to give rise to a greater 
molecular fraction. There are other effects working in the 
same direction. Larger grains are formed, and maintained 
unfragmented, preferentially in regions with recent massive 
star formation, and as the reaction from HI to H2 proceeds 
more rapidly under conditions of greater hydrostatic pres- 
sure (lElmegreenlll989l : iBlitz fc Rosolowskvl[2003 : iBlitz et all 
l2006h . and supernovae arising from massive stars give pe- 
riodic bursts of high pressure, these two effects combine to 
yield more massive molecular clouds in regions with previous 
high mass stellar populations. Following from this, if high 
mass molecular clouds yield stars with higher ranges of mass, 
we could expect to find a relation between dust grain size 
and the stellar mass range, because as we have mentioned 
above, lower mass stars tend to produce dust grains of lower 
size, which "seed" molecular clouds of lower mass, which in 
their turn produce a lower range of stellar masses. However 
there is a further effect which must be taken into account, 
which is that the grain distribution affects the internal ra- 
diative equilibrium of the molecular clouds and hence their 
tendency to condense into cores and then stars. The external 
radiating surfaces of the larger dust grains are proportion- 
ally smaller than those of their smaller counterparts with 
equal volume, so that larger clouds cool less efficiently and 
the Jeans masses of their cores are higher. This effect will 
tend to augment the tendency of the higher mass molecular 
clouds to yield higher ranges of stellar masses. The model we 
use therefore has two basic inputs: the cloud mass distribu- 
tion and the grain size distribution. In our first model (model 
A), we have taken a numerically derived MCGMF, from sec- 
tion 3.1.2 and folded it with the GSDF from section 3.1.1. 
In order to perform this folding we normalize the distribu- 
tions: we take the intervals in logr{fj,m) between -3.6 and 
-0.4 for the GSDF (see Fig. 2) and identify these extremes 
with those of the MCGMF: 4 and 6.5 in logMQ (see Fig. 
3) . Then we divide both intervals in 100 equal parts. So, we 
normalize dividing each linear value of GSDF by its higher 
value in the interval, and the same for MCGMF, then trans- 
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model A 

Kroupa (2001) 

PDMFKroupa(1993) 

Maciel & Rocha-Pinto (1998, variable SFR) 
Parravano et al. (2010) 
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Figure 4. Prediction of our model A, using a numerical molecu- 
lar cloud mass function (Section 3.1.2) folded with the dust grain 
size function (Section 3.1.1 , Fig. 2), compared with three obser- 
vationally based IMF's. 



the dust grains in the cloud, which is M^25r . This trend, 
at least in qualitative terms, can be understood from the 
known requireme nt of gra in surfaces for t he catalytic conver- 
sion o f HI to H2 jGould fc Sal peter 1 9631 : ICazaux fc Spaan3 



I2OO4I : IVidah et al.ll200d ). This implies that the mass of a 
cloud should increase with grain surface available. It is also 
known that bec a use o f the fractal geometry of dust grains 
iMutschke erall (|2009l ), who tested grain shape models in 
laboratory tests) the catalytic surface varies both with ra- 
dius and with shape. Spheroidal grains have effective sur- 
faces which increase as r^, so their surface per unit mass 
falls with increasing radius, obeying an r~^ relation. Toroidal 
shapes with constant section have effective surface and vol- 
ume varying as r^, so their effective surface area per unit 
mass is constant, and independent of radius. Fractal grains 
have effective surfaces which vary with their radii accord- 
ing to a power law with index very different from 2, which 
would lead us to anticipate a possible dependence of cloud 
mass on grain size with a positive index, as we can infer from 
the result of our normalization. It is worth noting here that 
the effective surface area for catalytic reactions will be, for 
fractal grains, much larger than the effective surface area for 
radiative equilibrium, since a major proportion of the frac- 
tal surface radiates into the grain, and only the outermost 
surface radiates energy away from the grain. Essentially the 
area to use in radiative equlibrium calculations will be ap- 
proximately that of a spherical grain with the same global 
radius, and will thus vary nearly as r^ even for a fractal 
grain. 



forming both distributions to other distributions normalized 
to unity in linear values. We then interpret the normalized 
MCGMF and GSDF as probability densities, multiplying 
them together point by point, and normalizing the results 
to 2 in log X(log M). Finally, we assimilate the extremes of 
both intervals (that of the GSDF and that of the MCGMF) 
to the extremes of a new interval of stellar masses at birth 
taken from 0.1 M© to 25 Mq. The curve we find using this 
semi-analytical approach, model A, is shown in Fig. 4 (full 
line). To explore an alternative approach we have used a the- 
oretical model for the MCGMF, based on CB02, in which 
the cloud mass distribution was computed using a simplified 
model in which the internal cloud dynamics was determined 
by the interaction of gravitational bounding and turbulent 
gas pressure. We fold this function with the theoretical grain 
size distribution of CBIO as in model A and the result is our 
model B (full line in Fig. 5) where we have extrapolated 
linearly the results to 100 Mq. 

In the process of identification of ranges (in fact the 
limiting values) of cloud mass and grain size we obtain two 
semi-empirical relationships. We then make a linear fit to 
a statistically representative set of 100 values of logr as a 
function of the 100 values of logMc and also for the 100 
values of the stellar masses at birth logM and so we obtain 
relations whose validity rests only on the good fits to data of 
the IMFs obtained from these relations (see Figs. 4 and 5). 
The first relation implies a functional dependence between 
the characteristic stellar mass, M, formed as a result of the 
internal collapse of a major gas cloud and the mass, Mc, 
of the cloud, and is just Mc^lO^M. The second is a rela- 
tion between the characteristic stellar mass, M, (in units of 
solar mass) and the characteristic linear size, r (in /im) of 



3.2 Analytical treatment 

Our model C is a theoretical model based on the assumption 
of balance between dust thermal emissi on and heating due to 
collisions with gas particles (Sch neider et al-l l|2006i) ). This 
leads to the equilibrium equation: 

AaT^Kpfj^scpgr = ngr2k{TG -Tgr)nHO gr ( ) /(lO) 

\TTmH / 

where a is the Stefan-Boltzmann constant, Tgr is the 
dust temperature, Ta is the gas temperature, Kp is the 
Planck mean opacity of dust grains, per unit mass, /3esc is 
the photon escape probability, pgr is the dust mass density, 
the grain number density, agr is the mean grain cross 
section, nn is the hydrogen number density, mn the mass 
of the hydrogen atom, and k is the Boltzman constant. The 
factor / takes into account the contribution to the gas-dust 
equilibrium of species other than hydrogen. If we take as 
an approximation that the gas temperature is significantly 
higher than the grain temperature, we can use Eq. (2) to 
derive a relatively simple relationship between the gas tem- 
perature and the grain parameters: 

Ty^cccjg^ng^ (11) 

As the mean grain cross section a^rocrg^, and the number 
of grains UgrXr'g'^ this gives us a relationship between gas 
temperature and grain size: 

Ta^^xrlr (12) 

As the Jeans mass within a cloud is proportional to 
t'c^Pc^^^ we find that the Jeans mass is proportional to Vg^. 
This result gives a clue to the possible causal link between 
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Figure 5. Prediction of our model B, which uses a theoretically 
derived molecular cloud mass function (Section 3.1.4) folded with 
the dust grain size function (Fig. 2), compared with the three 
chosen observational IMF's. 

the characteristic dust grain size and the final characteristic 
stellar mass, which we show gives a stellar mass- grain radius 
relationship where Mocr^''^ . 

We can now go on to obtain an analytic approxima- 
tion to the IMF, again assuming IMFcxGSDFxMCGMF, 
but now using an analytic fit to the data from CBIO, in 
the form GSDFcxr"^ '', and taking MCGMF as derived in 
section 3.1.3. This formulation gives us: 

IMFx]vr-^-''[e —+e ^ ] (13) 

If we now use, following section 3.1.3, »70cTg and starting 
from Eq. (12), and the relation between the Jeans mass and 
r, we find: 

IMFocM'^-^-'^[e ~+e ] (14) 

where a, rj and rj' are constants. 

4 COMPARISON OF THE MODEL 

PREDICTIONS WITH OBSERVATIONS 

We can see in Fig. 1 a comparison between the classical 
Salpeter's IMF and the PDMF of Kroupa ct al. (1 993) and 
also t he IMFs of: iKroupal lj200lh . iMaciel fc RochaTPintd 
ljl998( ) (in both cases, assuming co nstant SFR and variable 
SFR), and IParravano et all l|201ll '). In Figs. 4, 5 and 6 we 
can see the same data sample but compared with the pre- 
dictions of our models presented in this paper (A, B, and 
C respectively). Our models A and C can fit reasonably 
well all the data between 0.1 Mq and 25 Mq except that of 
IParravano et ID Iml) for masses lower than 1 Mq where 
the uncertainties in the observational data are largest. But 



Figure 6. Prediction of our model C, in which the grain size 
distribution is derived from a model where the gas temperature 
is maintained close to equilibrium by radiation from the grains, 
which thereby influence the Jeans mass within the cloud, and 
hence the stellar mass. The IMF is determined using analytical 
expressions containing the cloud mass and grain size distributions. 
The observational IMF's are plotted for comparison. 

the best fit is obtained for our model B where the GSDF is 
numerically derived and the MCGMF is derived analitically. 
We are not claiming here that we are in a position to decide 
finally which of the models gives the closest approach to the 
observations, particularly as the observations themselves are 
being continually revised. Nevertheless the introduction of 
the grain size dependence gives fits to the predicted IMF 
which are significant improvements on the use of the gas 
cloud mass function alone. We can see that the model A re- 
produces best the IMF derived by the later authors assuming 
a specific time-variable SFR produced from observational 
data of the chromospheric activity of local late type dwarfs. 
This agreement is in some sense model-dependent, because 
one of the underlying assumptions of the dust size distribu- 
tion component of the model was the same time variation 
sequence in the SFR. 



5 CONCLUSIONS 

This study offers prima facie evidence that there is a func- 
tional relationship between the characteristic masses of the 
stars at birth and the characteristic sizes of the dust grains 
which populate the molecular clouds giving rise to the stars. 
This relationship is superposed on a more conventionally 
accepted dependence of the stellar mass range on the pla- 
cental molecular cloud mass range. The evidence comes from 
the superior fits to the observationally derived IMF's of the 
models in which the two distributions, the MCGMF and 
the DGSF are folded together, compared to models based 
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on the MCGMF alone. However we have also offered a semi- 
quantitative explanation based on scenarios describing the 
effects of the dust grains on the formation of the molecular 
clouds, and on the collapse of the cloud cores to form stars. 
The details of these processes include two specific properties 
of the grains which might appear strange. In their catalytic 
action leading to the formation of molecular from atomic 
hydrogen, the formation rates favour larger grains, because 
their effective surface areas are fractal, so that their surface 
to volume ratio increases with grain radius. The same grains 
acting to radiate away heat show the opposite behaviour; 
their outer surfaces from which radiation can escape are 
more effective for the smaller grains, because the ratio of 
the radiative surface to grain volume falls with grain radius. 
These effects are both present in the intervention of grains 
in the cloud forming process and the eventual star forming 
process, and we have taken them quantitatively into account 
when deriving the IMS's in all the models presented here. 
The relation of dust grain size r to stellar characteristic mass 
M can be summarized in the expression M~25r°'^, where M 
is in units of solar masses, and r is in microns. The depen- 
dence of the characteristic stellar mass on the mass of the 
placental molecular cloud Mc can also be parametrized, and 
takes the form Mc^W^M. 
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